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Abstract 

We calculate the two-loop QCD correction to the form factors of on-shell b- 
quark decay to an energetic massless quark, which constitutes the last missing 
piece required for an 0{a1) determination of \Vub\ from inclusive semi-leptonic 
B — > Xuii^ decays in the shape-function region. 
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1 Introduction 



The strength of 6 — > m transitions, measured by the CKM element \Vub\, cahbrates 
one of the sides of the unitarity triangle and is therefore an important input to flavour 
physics. Exclusive semi-leptonic B —>■ Miu decays (with M = n,p, . . .) provide direct 
access to the b u transition, but require knowledge of the heavy-to-light meson form 
factors to extract \Vub\- Inclusive semi-leptonic heavy quark decays on the other hand, 
can be calculated in perturbation theory. However, the need to separate B —>■ Xuiu 
from an overwhelming charm final-state background requires cuts on the differential 
decay spectra, which in one way or another constrain the hadronic final state to have 
small invariant mass and large energy. The final state distribution then depends on a 
non-perturbative function describing the light-cone residual momentum distribution of 
the 6-quark in the B meson, the shape function [1,2]. Somewhat unfortunately, the 
two methods currently result in values of \Vub\ that seem to systematically differ, with 
the exclusive method favouring smaller \Vub\- This provides motivation and urgency to 
improving the theoretical accuracy of both methods. 

The differential distributions in inclusive semi-leptonic decay are calculated from the 
discontinuity 



W^" = - Im 

IT 



{B{pB)\t d'xe''^--T{j^^{Q)f{x))\B{pB)) 



(1) 



of the forward matrix element of a correlation function of the weak interaction current 
= u'y^{l — 75)^. In the shape-function region the structure functions into which the 
hadronic tensor W^'^ can be decomposed, can be organized into the factorized expres- 
sion [3] 

Wk = Yl ^ C.i^- ■ P) CM- -V) jdu J{pI) S{uj) (2) 

at leading power in the 1/mb expansion. The factors CiCj, J, S arise from the different 
scales nib, ^/"t^Aqcd and Aqcd, respectively, such that Ci and J can be computed in 
QCD perturbation theory, while the B meson shape function is non-perturbative. [a^ 
are numerical constants.) Eq. ([2]) has been worked out at order 0{as) in [4, 5] and 
forms the basis of the inclusive \Vub\ analysis performed in [6]. A summary of other 
inclusive \Vub\ analyses and recent experimental results based on these methods is given 
in [7]. At the two-loop order, the jet function J and the partonic shape function of the b 
quark have already been calculated [8,9]. Here we compute the two-loop hard matching 
coefficients Ci associated with the form factors of on- shell 6-quark decay into an energetic 
massless quark and provide a numerical estimate of the new term. All items for a 0{al) 
determination of \Vub\ from inclusive semi-leptonic decay are now in place, which should 
remove most of the perturbative theoretical uncertainty. A detailed phenomenological 
analysis is, however, beyond the scope of the current work. 

While this paper was in preparation, the calculation reported here has also been 
completed by Bonciani and Ferroglia [10]. The results have been compared prior to 
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publication and complete agreement has been found. On the day of submission of the 
present work, the paper [11] on the same topic has appeared. 

2 Structure of the calculation 
2.1 Set-up of the matching calculation 

The calculation of the short- distance coefficients Ci in ([2]) amounts to matching the 
current ^7^(1 — 75)6 to a set of leading-power currents in soft-coUinear effective theory 
(SCET) [12-14]. One-loop results for the have been obtained in this framework 
in [12,15]. The current is represented in SCET by 

[W7^(l - 75)&](0) = / dsY,C^{s) {iW,){sn+)V^,hM + . . . , (3) 

i 

where ^ is the coUinear up-quark field in SCET, the static heavy-quark field, and 
Wc a coUinear Wilson line. (We use notation as defined in more detail in [15].) There 
are three independent Dirac matrices that can appear on the right-hand side, which we 
choose as 

1^ = 7^^(1-75), r^ = t;'^(l+75), r^ = n^i(l + 75). (4) 

The ellipses in (|3]) denote higher-dimensional operators, which are not relevant to the 
present work, and s = smb. SCET makes extensive use of two light-like vectors n'^, 
with n+ ■ n_ = 2, with respect to which four-vectors are decomposed as 

= n+-p— + n^-p^+p1, (5) 

with n_ ■p± = n+ ■p± = 0. The coUinear field ^ describes modes which have n+ -p large, of 
order m^. To define the heavy-quark field, we choose a frame such that v'^ = (n^ + n^i)/2, 
that is n+ ■ f = n_ ■ w = 1. The factorization formula ([2]) for the differential decay 
distributions makes use of the momentum-space short- distance coefficient rather than 
the position-space expression appearing in the convolution in The momentum space 
coefficient functions are related to those defined above by 

a{u) = j dse'^'C,{s), (6) 

where the new variable u G [0, 1] equals the momentum fraction n+ -p/mb of the external 
up-quark line in a momentum-space Feynman diagram. 

The actual matching calculation is also done in momentum space and yields the 
momentum-space coefficient functions directly. To this end we consider the matrix ele- 
ment of the current between a bottom quark of mass rrtb and momentum pb and a massless 
up-quark with momentum p. Both quarks are considered on-shell, hence pi = ml and 
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= 0. The absence of any perturbative infrared scale implies that the SCET loop 
diagrams are scaleless, and the SCET matrix element is given by its tree-level expression 
multiplied by the universal renormalization factor for the SCET currents {^Wc)T'^h^'- 

{u{p) I i^W,) {sn+)T^K (OMPb)) = e^^"+-^ Zj m„_ Tf . (7) 

Thus, calculating the QCD matrix element in dimensional regularization yields directly 
the dimensionally regularized short-distance coefficient Ci{u). To this end, we first de- 
compose the matrix element of the hadronic current into three form factors according 
to 

{uip)\uY{l - I5)b\b{pb)) = F,{u) u{p)Y{l - 7rMPb) + ^2(w) u{p) — {^ + 'yMPb) 

nil) 

+Fs{u)u{p)^{l + j,)u{pb). (8) 
Pb-P 

The form factors Fi{u), i = 1,2,3 can only depend on the dimensionless variable u = 
2ph ■ p/ml and logarithms of /i^/m^, where /i is at once the renormalization scale of 
the strong coupling and the infrared factorization scale, since the on-shell heavy-to-light 
form factors contain soft and collinear divergences. Only Fi is non-zero at tree-level, 
Fi = 1 + 0{as), while ^2^3 = 0{as). Identifying p'^ = m^v'^ and = nij,p'^/{pb ■ p), 
we see that u equals the variable u defined in ([6]). At leading order in the heavy-quark 
expansion the spinors of the quark fields equal u{p) = u{pb) = u^, where the 

collinear and heavy quark spinors satisfy rjb-Un_ = and tfu^ = Uy, respectively. Eq. ([8]) 
then becomes 

3 

{u{p)\ur{l-l5)b\b{pb)) = J2Fi{u)un^T';uy, (9) 

1=1 

while using ([7]) the matrix element of the right-hand side of ([3]) equals 

3 

^ZjCi(?2+ -p/nib) (10) 

1=1 

so we simply have Ci{u) = Zj^Fi{u), i = 1, 2, 3. In the following we briefiy describe the 
method of calculating the two-loop QCD correction to the Fi. 

2.2 Calculational methods 

The computation of the two-loop QCD vertex corrections to semi-leptonic h ^ u decays 
involves the evaluation of the diagrams shown in Fig. [TJ We work in dimensional regu- 
larization with D = 4 — 2e, where UV and IR (soft and collinear) divergences appear as 
poles of up to the fourth order in e. For the matrix 75 we adopt the naive dimensional 
regularization (NDR) scheme with anticommuting 75. Furthermore, we treat all quarks 
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(b) (c) 




Figure 1: Two-loop diagrams needed for the calculation. Double straight lines stand for 
massive quarks of mass rub, whereas single ones stand for massless quarks. The filled 
circle in the last diagram represents the complete one- loop gluon self-energy with nf 
massless quarks (including the charm) and the massive bottom quark. 

except the bottom quark as massless. We checked that it is possible to include the effects 
of a non-zero charm mass analytically. However, since in the two-loop calculation of the 
jet and shape function [8,9] the charm mass is neglected, we also set it to zero in the 
present work. 

The amplitude of the diagrams is reduced by techniques that have become standard 
in multi-loop calculations. We apply a Passarino - Veltman [16] reduction of the vector 
and tensor integrals. The Dirac and color algebra is then performed by means of an in- 
house Mathematica routine. The dimensionally regularized scalar integrals are further 
reduced to a small set of master integrals (depicted in Fig. H] in the Appendix) using the 
Laporta algorithm [17,18] based on integration-by-parts (IBP) identities [19,20]. To this 
end we use the Maple package AIR [21]. 

The techniques we apply during the evaluation of the master integrals are manifold. 
The easier integrals can be written in a closed form in terms of F-functions and hyperge- 
ometric functions and subsequently expanded in e with the package HypExp [22,23]. In 
more complicated cases we derive Mellin-Barnes representations by means of the package 
AMBRE [24]. We perform the analytic continuation to e = with the package MB [25], 
which is also used for numerical cross-checks. We then apply Barnes' lemmas and the 
theorem of residues on the multiple Mellin-Barnes integrals, and insert integral represen- 
tations of hypergeometric functions as well as ■j/'-functions and Euler's i?-function where 
appropriate. As a third technique we apply the method of differential equations [26-28] 
and evaluate the boundary condition with the Mellin-Barnes technique. This renders, 
for instance, a three-dimensional MB representation in the case of the crossed six-line 
master integral of Fig. Hl^g) at m = 1. Eventually, the master integrals are evaluated as 
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Laurent series in e, with expansion coefficients of argument u expressed analytically in 
logarithms and polylogarithms of increasing weight. The maximum weight which ap- 
pears in our calculation is four, and we can express all functions but one in terms of 
ordinary polylogarithms. The function that cannot be expressed in terms of ordinary 
polylogarithms is the harmonic polylogarithm [29] HPL({— 2, 2}, l—u). We left, however, 
also a second function, HPL({— 1, 2}, l — u), in the HPL notation since its expression in 
terms of ordinary polylogarithms is rather complicated [30]. The master integrals have 
been calculated already in [31] and used in [32]. We find (almost) perfect agreement 
on the expressions in [31]. Moreover, we improve the numerical accuracy of one of the 
boundary conditions, see Appendix \K\ 



3 Renormalization 

After the pure two-loop calculation the result still contains divergences of UV as well 
as IR (soft and coUinear) nature. The former divergences are cancelled after addition of 
the UV counterterms, the latter disappear after inclusion of the jet and shape function 
contributions, since the partonic differential decay distributions are infrared-finite. In 
the following we also present the pole parts of the form factors Fi. Subtracting the IR 
poles as discussed in Sec. leads to the MS definition of the SCET current. 



3.1 UV renormalization 



In performing the UV renormalization we adopt the on-shell scheme for the heavy quark 
mass as well as for the heavy and light quark field. The strong coupling as on the 
other hand is renormalized in the MS scheme. The respective renormalization constants 
read [33-35] 



1-C, 



1-C, 



^ {D-l)T{2-f 



(47r)«/2 

13-4 



(D-S) 

(£)- i)r(2- 



(47r)«/2 



' (47r)^ ^ 2"^ 



Op 



51 



433 

+— - 137r2 + IGtt^ ln(2) - 24(3 + 0{p - 4) 



8 



(D-4)2 2(D-4) 

+ Cp Ca 



101 



803 



2{D-A) 
+CFnf tf 



+ 57r2 - Stt^ ln(2) + I2C3 + 0{D - 4) 
o 

8 18 47r2 59 



(11) 



22 



(D - Ay 



5 



+CFt 



f 



16 



38 



_(L>-4)2 3(L>-4) 



15^! + H5 + o(z)-4) 



(12) 



Here go is the bare QCD coupling. Cp = {N'^ — 1)/{2N) and Ca = N are the Casimir 
operators of the fundamental and adjoint representation of SU(A^), respectively, and 
tf = 1/2 denotes the normalization of the trace of two fundamental generators, nf 
stands for the number of massless quarks, and we set the number of heavy quarks of 
mass m to unity throughout the paper. The renormalization constant for a massless 
quark field in the on-shell scheme receives corrections only from two-loop and higher due 
to diagrams with massive quark loops. At two loops only a single diagram contributes 
and the Z factor assumes the following closed form 



1 + 2 t 



/ 



i)r(4-f)r(-f) 



{An 



(13) 



The renormalization constant of the strong coupling in the MS scheme reads 

a. II ^ 4 , ^ , 2 



1 + 



47r 



11 4 

-jCA--tf{nf + l] 



{D-A) 



(14) 



The renormalization of the quark fields amounts to mere multiplications. Expanding 
the bare and renormalized amplitude as well as the renormalization constants according 
to 



A-haxe — A^^'' + Qq A\^1^^ + Qq A\^1^^ + 0(5'o) 
-^ren = A^ + ^^en + fi'o '^ren + ^(S'o) ' 

Z, = l+gl5Zf^+gt5zf^ + 0{gl), 
we find for the renormalized amplitude 



i = m,h, I, 



A 



(1) 

bare 



A 



(2) 

bare 2 ^ bare 



2 



liszu^r+'-szi^^' 



(15) 
(16) 
(17) 

(18) 
(19) 



The bare coupling = Qq/IAtc) is then renormalized simply by the procedure a° = 
ZaOtslJ^^ with fj? = /i^exp(7£; — ln47r). This can also be seen from the way we present 
our results in ( l22l) . where Za accounts for the renormalization of the coupling constant. 

The only non-trivial contribution to the UV renormalization is therefore the one-loop 
diagram in Fig. [2l To this end, only that part of the counterterm Feynman rule that 
contains the one-loop correction 6Zm to Z^ has to be inserted, and the contribution 
has to be added to the RHS of Eq. (fT9l) . 
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Figure 2: One-loop counterterm diagram for mass renormalization 



3.2 IR subtraction, jet and shape function contribution 

Since the partonic structure functions Wi in ([2]) are infrared-finite, the 1/e poles of the 
two-loop coefficients Ci or, equivalently, Fi can be constructed independently using the 
two-loop expressions of the jet and shape functions [8, 9] and the one-loop coefficients 
Fi including their 0{e) parts (given below). To check our calculation, we perform the 
convolution J doj J{p1j) S{uj) of the unrenormalized jet and shape function with 0{al) 
accuracy and determine the infrared pole part of each Fi individually by forming appro- 
priate combinations of We find complete agreement with the poles obtained in the 
direct two- loop calculation of the Fj. 

In order to verify the pole cancellation in ([2]), one must remember that the calculation 
of the short-distance coefficients Ci is performed in a theory with five active flavours (the 
flfth being the bottom quark), hence as that appears above is ai^'*, while the jet and 
shape function are computed in four-flavour SCET. In combining Ci, J and 5* according 
to ([2]), the perturbative expansion of Ci must be expressed in terms of a^J^\ The D- 
dimensional relation between the renormalized couplings required for this purpose reads 
to one-loop accuracy 



a 



(5) 



^5 a 



(4) 



where ^45 = 1 + 



An 3 



Tie) 



nil 



(20) 



and the couplings are evaluated at the scale /i. 

The MS renormalized coefficient functions are obtained from Ci = Zj^Fi (see above), 
where [12] 



1 + 



1 ^ 1 



In 



(21) 



at one loop, and the two-loop renormalization factor can be obtained from the result of 
our calculation by requiring that Ci is free of 1/e poles at 0{al). 



4 Results 

We present the result for the short- distance coefficients (form factors) in the form 

F. = ) + ^ ) + ^ ) + Oia!) , (22) 
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with 



f(°) = 1, F^''> = 0, F^'" = 0. (23) 

The factor is given in f|T^ . and as well as in (122|) refer to a theory with five 
active quark flavors, Uf = A massless ones and one massive one of mass rrib- Note that we 
expressed the series in terms of ZaC^s rather than the renormalized coupling as- Thus, 
the contain 1/e poles that are cancelled by charge renormalization as well as the IR 
poles discussed above. 



.(0) 



.(0) 



4.1 One- loop results 

The one-loop hard matching coefficients are available in the literature through order 
O (e°) [5,12,15], but the two-loop analysis of the structure functions Wi requires them 
up to order 0{e^). Here we summarize the corresponding expressions for the form 
factors. 
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4.2 Two- loop results 

The two-loop part of the form factor Fi reads 
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21601^3 

4a2(,/ - l)lii(2) 57r2 



3 



u 



2(3 1327 

+ -77^ + 2co, 



A{u-1) u-1 16 



11, , , TT^ 58 



(42) 
(43) 



22, 2/ N 22, , , , , , 22 ^ , , llln(ii) 7r% 
-- ln\u) + - ln{u) ln(l - u) + - U,{u) + ^^j^ - y H'^) 

166, , , 11, 13l7r2 6301 
+-9-ln(H) + -C3-^-^, 



(44) 



(4^2 - lOu + 7) 

3^l3 



+24Li4(ii) + 24Li4 



\n^{l-u) + ln2(l - m) - 41n^(l - it) ln(H) 
2 {7u^ - 16^2 + 40m - 28) 



u 

u - 1 



ln(ii)C3 



(4^3 + 12^2 - 30m + 21) , 2/ M 2/1 N 

^ ^ W{u) ln2(l - w) 



(106^3 + 33^2 - 84m - 9) 2 



27r2 



(4^2 + 4^2 - 10m + 7) 



ln2(M) ln(l - m) + — ln2(M) + 



[Li2(M) + 2 ln(M) ln(l - u) \a2{u)\ 



221n(l -m) ln(M) 
3(m- 1) 



{w- + 2m - 2) 



27r2 



8 HPL({-2, 2}, 1 - m) - — Li2(^/ - 1) 



+ 



171n2(M) 



4 (5x^2 - 10m + 4) 



TT 



HPL({-l,2},l-M) + -ln(2-M) 

4 TTT.T/r . . N 7r2ln(2-M) 88 , , 
— HPL({-1, 2}, 1 - «) + ^ji-^i + - 



2u 

547 
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ln"(M) 



8a2(1,/2 - ll» + 8) , , , , 4129, , , 351n2(M) 

+ - Hu) + ln(«) - ^^^^ 



(215x^2 - 108ii+ 153) 



9^/2 



[ln(l - m) ln(M) + U2{u)] 



27r2ln(M) 5331n(M) 
S{u-1) ^ lS{u-l) 



{u^ - 3) 



\[t?{u) ln(l + -u) + TT^ ln(l + -u) + 2 ln(M) U2{-u) - 2 Li3(-M)] 

3 (2m3 + 8^2 - 14m - 1) 



2m3 

7r2(155ii2 - 180ii + 216) 



18x^2 



ln(M) 



■ \n.{u) Li2(tt) 
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+ + 3^^3 Li.(«) Hu) Li3(l - u) 

(2u^ + 93u'^ - 90u + 9) ^ , , (70ii=^ - lllii^ + 168ii + 9) , 

Li3(l -u) —. Li3(m) 



2 (20x^2 - 58ii + 43) ^ , , tt^ (230x^2 - 405ii + 378) 
TT^ (269it=^ + 4801*2 - 1236w + 876) ^ (461w=^ - 396it2 + 1008it + 54)C3 



10801*3 181*3 
2 7r2(it-4)ln(2) IItt^ 146461 

20 4 

P-sH^y-glnH, (46) 

, , 8, 2/ ^ 8, , ,^ , 8_^. , , 41n(ii) 56, , , ISvr^ 557 
p_.(.) = - X.\u) - - Hu) ln(l - .) - - U,iu) - 3^ - - ln(.) + ^ + ^ , 

(47) 

32, 3, , 16, 2, N 81n=^(M) 112 W(m) 8 ln(M) ln(l - m) 



9 ' ' 3 ' ' ' ' 3(xi-l) 9 3(i/-l) 

112, , , , 861n(M) 267r2, ^ ^ 706 , 8Li2(M) 74, 
ln{u) ln(l - -u) - — — — ln{u) - — ln{u) - ^ . ^ ' - —(3 



9 ' ' ' ' 9(m-1) 9 ' ' 27 ' ' 3(m-1) 9 

112^.,, 16,.,, , 32,.,, 47r2 1067r2 11813 

Li2(w) + — Li3(l -u) + — Li3(m) + — — + — — + — — , (48) 



9 ' 3 ' 3 ' 9(m- 1) 27 324 

M^^) = y-^M^), (49) 

g_i(K) = ^ ln2(^) - ^ In(^) ln(l - t.) - ^ U,{u) - ^^^^ - 4 Hu) + ^ + 8 , (50) 

/ \ 16, 3, , 8, 2/ M /-, X AH^iu) 2/ \ 81n(M) ln(l — m) 
,„(«) = -- ln»(„) + - ln^«) 1„(1 -u) + + 4 ln^(„) ^'^^^ 

8 (Sti^ - 39u2 + 54u - 16) ,, , , , ,, , , ,, 861n(ii) 4ir\ , ^ 

2 (409^2- 660^ + 192), , , 8Li2(M) 16(^3-3^ + 3), , 
^ Mn(ii)-— — ^ + — ^ — — ^Li3(l-'u) 



27^2 ^ ^ 3(m-1) 3^3 

16, , , 47r2 7r2(llM2 + 12m - 24) 16 (5'u3 - 9^ + 9) , 

+T ^''(''' + - -^^^ <^ 
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We comment on the role of the constant cq in the Appendix. The two-loop part of the 
form factor F2 reads 

^ = Cl 1^ r_2(K) + ^ [2Lr-2{u) + r-i{u)] + [21" r_2(K) + 2Lr_i(K) + ro{u)\ | 



+Cf Ca \ - s-i{u) + [2L s-i{u) + soN] 



+CFtfnf l^i-iH + [2Li_i(ii) +toH]| 



with 



r-2{u) = -go{u) , 
2u 

ro{u) 



{u - 1)2 
S (2m - 7) 



[Li.(l-.) + 31nV)]-iH^^4^+ ' 



U-l 

ln^(l - li) + 277^ ln2(l - li) - 41n^(l - li) In(ii) 



(52) 

(53) 
(54) 



-24 \vl{u) C3 + 9 In^(w) ln^(l - -u) + 3 Li|(ii) + 6 hi{u) ln(l - -u) U2{u) 



+24 ln(M) Li3(l - -u) + 24 Ui{u) + 24 Li4 



2 (11^2 + 52m -3) 



ln^(ii)ln(l -m) + 



27r2 



M- 1 

(23m - 36) 



+ 



50W(m) 



M 



[ln(l - m) ln(M) + Li2(M)] 



8 HPL({-2, 2}, 1 - m) - — Li2(M - 1) 



2(3m-8), 2t ^^ m ^ 

-A— ^inV)m-x^) 



32 {u + 2) 

V? 
16 M 



TT 



HPL({-l,2},l-M) + -ln(2-M) 



(25^2 - 49m + 28) , 2 



(m-1)^ 



In'(M) 



(m - 1)2 

327r2 (4m - 13) 
3^l3 



HPL({-1, 2}, 1 - m) + ^ ln(2 - u) + ^ \n\u) 



ln(l - m) ln(M) - ^^^^""^ [ln(l - m) ln(M) + Li2(M)] 



2(M2+4M+3)r, 2 



[in^(M) ln(l + m) + TT^ ln(l + m) +2 ln(M) Li2(-M) - 2Li3(-M)] 
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, 327r^(i^ + 6) , , , , A(17u' + 80u-3) ^ , ^TW ^ 12(3x^-4) , 
H TT^ ln(M) H In(ti) Li2(ti) In(ti) Li2(ti) 

u"^ [u — ly 

^. , , 4(8m + 3) ^. ,^ , 2(9m-8) ^. ,^ 
Ll2 M + ^ 3 ' L13 1 - M - L13 1 - u) 

^. . . 4(15«-16) ^. , , 16(10m-31) ^. ,^ 
Usiu) + ^7 Lisiu) ^—^^ ^ Li4(l - u) 



3x^2 
647r2 (m - 3) 



3u^ 

4 (23x^2 + 108ii - 3) 



27r2 (25'u - 84) 71^ (41m - 54) 2%^ (40m - 137) 7r2(3lM - 32) 
H ttt; I r;^ 5 ;n ttt; lii(Mj 



3m2 

8 (3m - 1) 

~ [u-iy 

22u \n{u) 
3(m- 1)2 ~ 3(m- 1) 

2 (4m - 7) 



6(m- 1)2 



45 



3(m- 1)2 



, 4 (22^2 + 116m - 3)C3 4 (15m -16), 31 



M"" 



u-l 



22 



(55) 
(56) 



3m3 



ln^(l -m) + 27r=^ln^(l -m) -41n^(l -M)ln(M) 



-241n(M)C3 + 91n^(M)ln2(l -m) + 3Li2(M) + 61n(M) ln(l - m) Li2(M 



+24 ln(M) Li3(l - m) + 24 Li4(M) + 24 U4 



u 



u- 1 



+ 



17W(m) 



M 



(m2 - 32m - 3) 



M"' 



\u^{u) ln(l - m) + 



27r2 



4 (17m -6) 
3(m- 1)2 



[ln(l - m) ln(M) + Li2(M)] 



8 HPL({-2, 2}, 1 - m) - — Li2(M - 1) 



M 

16 (m + 2) 

m2 

8m 



M 



TT 



HPL({-l,2},l-M) + -ln(2-M) 



— ln2(M)ln(l-M) 



5 (13m -6), . 



3(m- 1)2 



(m - 1)2 

647r2 (m - 2) 
3^l3 

(m2 + 4m + 3) 



HPL({-l,2},l-M) + -ln(2-M) 



ln(l - m) ln(M) - ^^^^""^ -^^^ [ln(l - m) ln(M) + Li2(M)] 



[\n^{u) ln(l + m) + TT^ ln(l + m) + 2 ln(M) Li2(-M) - 2Li3(-M)] 



47r2(M+18) 
3^l2 



ln(M) 



2(2m2 - 46m - 3) 



\n{u) Li2(M) 



12 



ln(M) Li2(M) 



87r2 (4m - 9) , , , 2 (5m2 + 34m - 3) , , 
^— Li2(m) + ^ ^ ^ Li3(l - m) - 



3m=^ 



M"" 



M- 1 

4Li3(l-M) 
M- 1 
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6(u'^ -20u-l) ^ . , , 16 Ligfw) 4 (20ii - 43) , . 
f ^ ^ ^ Ll3 M + ^ - ^ ^ ^ Ll4 1 - u 



27r2 {?>u - 14) _ 27r2 (8ii + 3) tt^ {^u - 73) _ 47r2(3ii - 2) 
9(ii-l)2 ^ 45il3 3(ii- 1)2 



■ ln(M) 



(323m - 54) 2 (2x^2 - 52ii - 3)C3 _ 16 Cs _ 401 



8m ln(Ml 



'-^^^) = -3F^ + ^(^' 
^'^""^ ^ 9(^^-1)2 [l21n^W + 12Li2(l - m) - 191n(M)] + , 

8m n 2/ n ^t. M 4(79m- 60), , , 801n(M) 

[2 ln(M) ln(l - M) + 2 Li.(M) + n^] + 

^32Li3(l -m) _ 32C3 _ 104 

3u 

Finally, the two-loop part of the form factor reads 

=C||^m;_2(m) + ^[2Lm;_2(m)+m;-i(m)]+ [2L2 m;_2(m) + 2Lm;_i(m) + 

x_i{u) + Xo{u)] > 

+CFtf Uf |^|/_i(m) + [2Ly_i{u) +2/o(m)] 

+CFtf |^;2_i(m) + [2L z-i{u) + Zo{u)] 
with 

w-2{u) = -ho{u) , 

, , m(2m — 1) r , 9/ V ^ / vi m(24m — 17) , , , 9m 



4 (2m2 - 16m + 21) 



ln^(l - m) + 27r2 ln^(l -m) -41n-^(l -M)ln(M) 
-24 ln(M) Cs + 9 ln2(M) ln2(l - m) + 3 ^^(m) + 6 ln(M) ln(l - m) Li2(M 
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+241n(M) Li3(l -u) + 24Li4(M) + 24Li4 



u 



u-1 



u 



2{u-l) 



^ln^('u)ln(l-'u) + y^ -|ln('u) Li2(m) 



6(m-2) 



- 1)^ 



+:r7 TV^T [ln(l - u) \vl{u) + Li2(w)] - 7 \v?{u) ln(l - + 



2(m- 1)2 
32(2m - 3) 

,2 



(t.-l)^ 



TT 



HPL({-2,2},l-ii)--Li2(ii-l) 



16 {v? -u + 6) 



8 (3m - 2) 
(m-1)2 



TT 



HPL({-l,2},l-ii) + -ln(2-ii) 



2 (13m - 14) 
(m - 1)2 

7r2(29M-30) 



+ 



Li3(M) 
\vl{u) 



HPL({-1, 2}, 1 - m) + ^ ln(2 - «) + ^ ln'(«) 



6(m- 1)2 
+ Stt^ ln(2) 



(24^3 - 73^2 + 69m - 24) 



levr^ (4m2 - 30m + 39) 
3^l3 



ln(l — m) ln(M) 



n /1 NT / A T - / M 7r2(M2-80M + 288) , , , 
[ln(l - m) ln(M) + Li2(m)] ^ ln(M) 



2(m- 1)3 
_2(6m2 -71M + 75) 

m2 

(2m3 + 3m2 + 8m + 9) r, 2/ X, / ^ 9, / X , / ^. / X 

+^ [in^(M) ln(l + m) + TT^ ln(l + m) + 2 ln(M) Li2(-M) 

. ,^ 2(10m3-69m2 + 244m-9), , , ^. , , (11m -10)^.,^ 

-2L13 -M - ^ In M L12 M + ^- L13 1 - M 

M-^ [u — ly 

327r2 (xi2 _ 7^ + 9) 2 (3m3 + 8m2 - 28m - 9) , . 

3^l3 Li2(m) + Ll3(l - m) 

2 (4m3 - 91m2 + 328m - 9) 



M*" 



8(10m^-72m + 93) ^. 

Li3(m) ^ — Ll4(l - m) 



M*" 



TT^ (5m2 - 390m + 504) tt^ (17m - 30) 2 (14m3 - 94m2 + 344m - 9) (3 



6m2 



12(m- 1)^ 



M"" 



(79m2 - 67m + 4) , , , 7r^(40M2- 314m + 411) 2 (13m - 14) ^ 
+^^77 TTT^ ln(M) + ^ —-^ + -^Cs , (65) 



2(m- 1)2 
11m(2m- 1) 



45 m3 



(m - 1)^ 



ln(M) + 



11m 



3(m- 1)2 ' 3(m- 1) ' 

Iv? - 24m + 21) 



(66) 



3m3 



ln^(l - m) + 27r'ln'(l - m) - 41n3(l - m) ln(M) 



-24 ln(M) Ca + 9 In^(M) ln^(l - m) + 3 Li^(M) + 6 ln(M) ln(l - m) Li2(M) 
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y-i{u) 
yo(u) 

z-i{u) 



+241n(M) Li3(l -u) + 24Li4(M) + 24Li4 



u 



u-1 



51W(m) 347x^ + 54 
2^ ^ lS{u-l) 



{2u^ - 51u^ + 92u + 9) 
2{39u - 28) 



16{2u - 3) 



ln^{u) ln(l - m) H ln{u) U2{u) 

XL i 

[ln(l - u) hi{u) + U2{u)] + ln2(M) ln(l - u) 



u-l 



TT 



HPL({-2,2},l-u)--Li2(K-l) 



Li3(ii) 



8 (m^ - -u + 6) 



4 (3m - 2) 



TT 



HPL({-l,2},l-ii) + -ln(2-ti) 



u-l 

27r2(5?7, - 4) 



+ 



\n{u) 



TT 



HPL({-l,2},l-ii) + -ln(2-ii) 



(94^2 _ 53m - 6) 



16n^ (4^2 - I3u + 12) 



- An^ ln(2) 



\n{l-u) \n{u) 



6 {u - 1)2 

_ (56^2 -216m + 153) 
3^l2 

(2^3 + 3^2 + 8m + 9) 
2^l3 

. , (8^3 - 78^2 + 134m + 9) , , , , 2Li3(l-M) 
-2Li3(-M)] - ^ ^ ln(M) Li2(M) + ^ 



rwi ^^(^^J■(^^ 27r2(2M2 - 29m + 54) ^ ^ ^ 

ln(l - m) ln(M) + Li2(m) ^ — \n{u) 

3m^ 

[ln2(M) ln(l + m) + 7r2ln(l + m) + 21n(M) Li2(-M) 



M"" 



M- 1 



47r2 (8m2 - 28m + 27) ^ , , (6u^ - Alu^ + 98m - 9) , , 

Li2(m) Ll3(l - m) 



3m3 

(14m3 - 105m2 + 176m + 9) 



M"^ 



^ . , , 2 (40m2 - 136m + 129) ^ . ,^ , 
Li3(m) ^ Ll4(l - m) 



M"" 



7r2 (76m2 - 285m + 378) tt^ (30m - 19) {Au^ - 102m2 + 160m + 9)C3 



9m2 



9(m- 1)^ 



M-" 



(760m2 - 329m - 162) ^^^^^ ^ tt^ (80m2 - 246m + 219) ^ 8(3 



18 (m- 1)2 
4m(2m - 1) 



90 M'^ 



M-1 ' 



3(m- 1)2 

_8m(2m- 1) 
~ 3(m- 1)2 

y-i{u) , 



\n(u) 



4m 



3(m-1) ' 



(67) 
(68) 



r, 2/ X T. M 2m(68m-49), , , 62m 

[ln2(M) + Li2(l - m)] + .^2 Hu) - Tw:— TT ' (69) 



9(m- 1)2 



9(m- 1) 



(70) 
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4ii(2u- 1) , J, ^ 8(3ii-2),.„ ^ 2(14711- 128),,, 



9u 

2(55m- 24) 47r2(2M2 -3m+ 18) 16(2m - 3) ^ , , 284 

H ^ Li3(l-M) + — . (71) 



9(m-1) 9m2 m3 / 3^ 

As already mentioned in the introduction, our results have been compared analytically 
with those of [10] and complete agreement has been obtained. 



5 Numerical evaluation and conclusion 

In Fig. [3] we show the coefficient functions Ci{u) in the one- (dashed) and two- loop 
(solid) approximation. For this purpose, we define the renormalized matching coefficients 
in the MS scheme by subtracting minimally the infrared poles. More precisely, from 
Ci = Zj^Fi, we obtain 



2 



a = G.'»' + ^cW + i^cf +0(a^). (72) 



with 



Ci'^ = F!'^ + [ZjTF!^. (74) 
Cf) = Fi') + ^(1) I [^-1] (1) + 4) + ^(1) } + [^-1] m ^(0) . (75) 

Expressions for the F^^"^ have been given in earlier sections. [Zj^]*^^) can be derived by 
inverting (!2T|) and taking the coefficient of of' /{An). ^I^-* is the coefficient of a'f' /{An) 
m ([20]), and zi^^ the one of Us/ {An) in f|T4l) . [Zj^Y^^ is defined to subtract the remaining 
pole parts in ^ + f[^\ [Zf]^^'^ + + zi^^}, such that Cf^ is IR-finite. Note that to 
evaluate fl75l) one needs the 0{e'^) terms of F^^"* given in Sec. 14. H since [Zj^]^^) contains 
a 1/e^ singularity. 

These coefficients depend on the IR factorization scale, which cancels only in the 
product of hard, jet- and shape-function factors. In physical applications the factoriza- 
tion scale ranges between ^/rr^AQ^ and mi,. To illustrate the size of the new two-loop 
correction, we therefore evaluate Ci{u) at /i = = 4.8 GeV (blue/dark grey curves) 
and at /i = 1.5 GeV (orange/light-grey curves). We emphasize that the difference be- 
tween these two choices is not a theoretical error - it is compensated by a correspond- 
ing scale dependence of the convolution J -k S. The value of the strong coupling is 
af\mh) = ai'^\mb) = 0.22, and ai''^(1.5 GeV) = 0.3753. Comparison of the dashed 
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Figure 3: The matching coefficients Ci{u) {i = 1,2,3 from top to bottom) in the one- 
loop (dashed) and two- loop (solid) approximation. The blue/dark grey curves refer to 
= mfe = 4.8 GeV, the orange/light grey curves to n = 1.5 GeV. 
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and solid curves of the same colour in Fig. [3] shows that the two-loop corrections are 
generally very moderate, if not small, except in the region of small u, where increasing 
powers of large logarithms take over. The implications of this result for the \Vub\ de- 
termination remain to be investigated. The impact of the 0{al) terms depends on the 
combination CiCj J-k S, and the numerical size of the two-loop correction to the jet and 
shape function has not yet been analyzed. A straightforward evaluation of the partonic 
structure functions Wi in the shape-function region indicates sizeable two-loop effects. A 
reanalysis of existing B XuCv decay distribution data with 0{a1) accuracy taking into 
account renormalization group summation and a model of the B meson shape function 
is therefore well motivated. 
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A Master integrals 

The two-loop calculation gives rise to 18 master integrals, which are depicted in Fig. |H 
All necessary master integrals were already computed in [31], which chooses a slightly 
different basis compared to the present work. Certain individual master integrals can 
also be found in [36-40]. We find almost perfect agreement with the results in [31]. In 
the case of the crossed six-line master integral Ig in Fig. 111(g) we improve the numerical 
accuracy of the boundary condition of the finite part. We therefore give our results for 
Ig explicitly. Our integration measure reads 




(76) 



We define the prefactor 



1 



(77) 



(47r) 



'/^r(i-e) 



and express our results in terms of 



X = 



{Pb - pf + if] 
ml 



l-u. 



(78) 
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Figure 4: Two-loop master integrals needed for the calculation. Double lines are massive, 
and single lines are massless. All diagrams stand for scalar integrals with unit numerator. 
Dots on lines represent squared propagators. Topologies with one or more dots stand for 
the undotted diagram and all diagrams with one single dot, i.e. topology (h) stands for 
three diagrams. 



Here +iri stems from the +iri prescription which we tacitly assume to be included in the 
propagators of the integral below. We find 



[(^2 + VbY - ml] {k2 + pf [{ki + - ml] (fcs -ki+ 
1 

X 



{k^-kifkl 



(-4) (-3) (-2) (-1) 



with 

1 



c 



(-4) 
9 



12il-x) 



22 



.(-3) 



.(-2) 



.(-1) 



.(0) 



ln(l — x) 
3(1 ' 

1 



72 (1 - x)2 
1 

36(1 -x)2 
1 



267 Ca] 



l-x) 



[48111^(1 -x) - Stt^] , 

[-321n^(l -x) + lOvr^ ln(l - 

8 5 

-ln^(l -x) - -vr^ 111^(1 -x) + 8 ln(l -x) Li3(x) 

9 9 

65 

+y ln(l -x)C3 + 4Li2(x) +co 



(80) 



The constant Cq, obtained with the package MB.m [25] from a three-dimensional Melhn- 
Barnes representation, equals —60.2493267(10), where the number in parenthesis gives 
the uncertainty of the last two digits displayed. The number excludes Cq = — 897r'^/144 
which was found in [3lJ^. It suggests Cq = — 1677r'^/270, in agreement with the proposal 
in [10]. To date there does, however, not exist a value for cq that is derived completely an- 
alytically. Except for cq, we obtained all other terms in the master integrals by analytical 
steps, that is, without fitting rational numbers to numerical values. 
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